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Abstract 



In the recent years we have shown that cylindrical biological membranes such as nerve 
^D axons under physiological conditions are able to support stable electromechanical 

pulses called solitons. These pulses share many similarities with the nervous impulse, 

^ e.g., the propagation velocity as well as the measured reversible heat production 

h^ and changes in thickness and length that cannot be explained with traditional nerve 

^H models. A necessary condition for solitary pulse propagation is the simultaneous 

\^ existence of nonlinearity and dispersion, i.e., the dependence of the speed of sound on 

density and frequency. A prerequisite for the nonlinearity is the presence of a chain 

melting transition close to physiological temperatures. The transition causes a density 
rS dependence of the elastic constants which can easily be determined by experiment. 

The frequency dependence is more difficult to determine. The typical time scale of a 

nerve pulse is 1 ms, corresponding to a characteristic frequency in the range up to one 

r^ kHz. Dispersion in the sub-kHz regime is difficult to measure due to the very long 

c/2 wave lengths involved. In this contribution we address theoretically the dispersion 

of the speed of sound in lipid membranes and relate it to experimentally accessible 



CZ3 relaxation times by using linear response theory. This ultimately leads to an extension 



of the differential equation for soliton propagation. 
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1 Introduction 

Biological membranes are ubiquitous in the living world. Despite their diversity in com- 
position, membranes of different cells or organelles are remarkably similar in structure 
and exhibit similar thermodynamic properties. They exist as thin, almost two-dimensional 
lipid bilayers whose primary function is to separate the interior of cells and organelles 
(sub-cellular compartment) from their external environments. This separation leads in turn 
to the creation of chemical and biological gradients which play a pivotal role in many 
cellular and sub-cellular processes, e.g. Adenosine Tri-Phosphate (ATP) production. A 
particularly important feature of biomembranes is the propagation of voltage signals in the 
axons of neurons which allows cells to communicate quickly over long distances, an ability 
that is vital for higher lifeforms such as animals U_, 2J. 

Biological membranes exhibit a phase-transition between an ordered and a disordered 
lipid phase near physiological conditions 1 3 1 . It has been shown that organisms alter their 
detailed lipid composition in order to maintain the temperature of this phase-transition 
despite different growth conditions |S^6l. The biological implications of membrane phase- 
transitions continue to be an area of active research. Near a phase transition the behavior of 
the membrane changes quite drastically: The thermodynamic susceptibilities, such as heat 
capacity and compressibility, display a maximum and the characteristic relaxation times of 
the membrane show a drastic slowing down IIt UTTI . 



The melting transition in lipid membrane is accompanied by a significant change of 
the lateral density by about -20%. Thus, the elastic constants are not only temperature 
dependent but they are also sensitive functions of density. Together with the observed 
frequency dependence of the elastic constants (dispersion) this leads to the possibility 
of localized solitary pulse (or soliton) propagation in biomembrane cylinders such as 
nerve axons. With the emergence of the Soliton theory for nerve pulse propagation, the 



investigation of sound propagation in lipid membranes close to the lipid melting transi- 
tion has become an important issue [2J. The Soliton model describes nerve signals as 
the propagation of adiabatic localized density pulses in the nerve axon membrane. This 
view is based on macroscopic thermodynamics arguments in contrast to the well-known 
Hodgkin-Huxley model for the action potential that is based on the non-adiabatic electrical 
properties of single protein molecules (ion channels). Using this alternative model, we have 
been able to make correct predictions regarding the propagation velocity of the nerve signal 
in myelinated nerves, along with a number of new predictions regarding the excitation of 
nerves and the role of general anesthetics lfT2l . In addition, the Soliton model explains 
a number of observations about nerve signal propagation which are not included in the 
Hodgkin and Huxley model, such as changes in the thickness of the membrane, changes in 
the length of the nerve and the existence of phase transition phenomena |[T3l . The solitary 
wave is a sound phenomenon which can take place in media displaying dispersion and 
non-linearity in the density. Both of these criteria are met close to the main lipid transition. 
However, the magnitude of dispersion in the frequency regime of interest for nerve pulses 
(up to 1 kHz) is unknown fH. Exploring sound propagation in lipid membranes is thus 
an important task for improving our understanding of mechanical pulse propagation in 
nerves. All previous attempts to explore sound propagation in lipid membranes have 
focused on the ultrasonic regime ||9j [T4l - [T6l . and it has clearly been demonstrated the 
dispersion exists in this frequency regime. Furthermore, the low frequency limit of the 
adiabatic compressibility of membranes (which determines the sound velocity) is equal to 
the isothermal compressibility, which is significantly larger than the compressibility in the 
Mhz regime. With the additional knowledge that relaxation times in biomembranes are of 
the order of milliseconds to seconds, it is quite plausible to expect significant dispersion 
effects in the frequent regime up to 1 kHz. 

Theoretical efforts to describe sound propagation in lipid membranes near the lipid melt- 
ing transition in the ultra-sonic regime have been based on scaling theory, which assumes 
critical relaxation behavior during the transition llT6in7l . However, a number low frequency 
experiments, pressure jump experiments lITOl [TSl and stationary perturbation techniques 
ifTTl [T9l all show non-critical relaxation dynamics. These findings have led us to propose 
a non-critical thermodynamical description of sound propagation in lipid membranes near 
the lipid melting transition for low frequencies based on linear response theory. 

In this article we will present a theoretical derivation of the magnitude of dispersion for 
membranes close to lipid melting transitions. The goal is to modify the wave equation for 
solitons in biomembranes. This will ultimately lead to a natural time scale for the pulse 
length, which we will explore in future work. 

2 The propagating soliton in nerve membranes 

In the following we present the hydrodynamic equations that govern the propagation of 
density waves in cylindrical membranes in general, and in nerve membranes close to the 



chain melting transition in particular. 

In it simplest formulation the wave equation for compressible fluids assumes the forrrQ 



where 



^'^ - V(c2Vp) (1) 




(2) 



5,0 ^/^ 

is the speed of sound for low amplitude waves (Ap <^ po)> ^5 is the adiabatic compress- 
ibility, and p{x, t) is the density. If the speed of sound is roughly independent of density, 
this equation simplifies to 

The wave equation in one dimension is then given by 

' c^^p . (4) 



dt"^ dx \ dx ^ 
For low amplitude sound we further assume that there is dispersion of the form 

c2 = cg + /iow2 + ..., (5) 

which corresponds to a Taylor expansion of the sound velocity with respect to frequency. 
The parameter /iq indicates the magnitude of the dispersion. Due to symmetry arguments, 
only even power terms appear in this expansion. One way to generate this frequency depen- 
dence is to add a dispersion term to the wave equation 

^ = i_fc2A \-h— . (6) 

dt"^ dx \ dx J dx'^ 

The density of a small amplitude plane wave can be written as 

p{x, t) = pq + Ap with A/9 = A sm(kx — cot) = A sm{k{x — ct)) . (7) 

The amplitude of this plane wave is A, and its velocity is with c = uj/k. Inserting this 
into Eq. ffl yields the dispersion relation in Eq. pj with /iq = h/cQ. We have shown 
experimentally that the sound velocity close to melting transitions in lipid membranes is a 
sensitive nonlinear function of density. Thus, we expand 

c^ = cl+pAp + q{Apf + .... (8) 

'a derivation of tlie equation of sound, based on fluid dynamics, can be found in [201 . There are two basic 
assumptions in the derivation of the equation of sound: Perturbations are small, and sound propagation is an 
adiabatic process. 



The parameters p and q describe the nonhnear elastic properties of membranes. At temper- 
atures slightly above the melting transition, lipid membranes have negative values for the 
parameter p and positive values for the parameter q. The final wave equation is given by 



5t2 



l(icl+pAp + qiApf)lp 






(9) 



We have shown that this equation possesses analytical solitary solutions that in many 
aspects resemble the nerve pulse (see Fig. [TJ. 
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Figure 1 : The propagating soliton using parameters appropriate for unilamellar DPPC vesicles 
and a dispersion constant h = 2 m^/s^ (from |21 1). The soliton has a width of about 10 cm and 
a duration of about 1 ms, which is very similar to action potentials in myeUnated nerves. 



While the above equation makes use of the fact that the speed of sound is a known 
function of density, the dispersion constant h must be regarded as an adjustable parameter 
due to the absence of quantitative empirical data regarding dispersion in the low frequency 
regime. The magnitude of h sets the width and the time scale of the mechanical pulse. In 
previous publications it was adjusted to /i = 2m^ / s^ in order to match the observed width 
of the nerve pulse, which is about 10 cm. However, we will argue below that h is expected to 
be density dependent and that its functional form can be approximated using experimental 
knowledge about relaxation time scales and elastic constants. This will ultimately lead to a 
wave equation for the mechanical pulse in nerve axons that is free of adjustable parameters 
and has a time scale that is fixed by the thermodynamics of the system. 



3 Brief Overview of Sound 

Sound is a propagating low amplitude density wave in compressible medium which, due to 
its adiabatic nature, is accompanied by a corresponding temperature wave. The equation 
governing sound propagation is universal. This generality implies that sound propagation is 
determined solely by the macroscopic thermodynamical properties of the system. 
As mentioned above, the equation of sound for low amplitude waves has the form: 



The general solution has the form: 



5^2 



Aexp(ztj(t- j;/c)), (10) 



which is merely eq. Q in complex notation. Due to dispersion and the absorption of sound 
in a real medium, the effective speed of sound, c, is a complex quantity. The real part of 
the speed of sound will have a phase shift (as a result of dispersion), the imaginary part will 
lead to a decrease in the amplitude or intensity of the sound as it propagates (attenuation). 
This can be seen by inserting the complex speed of sound into Eq. ([TO]). 



/9 = ylexp (i(^ (t — 2;i?e(c)/|c| )) exp (— xa;/?7i(c)/|c| ) , (11) 

where 



Reic) 
is the effective speed of sound which would be measured in experiment. 



(12) 



In 1928, Herzfeld and Rice extended the theory of sound by arguing that internal vi- 
brational modes of polyatomic molecules require time to approach thermal equilibrium 
with translational degrees of freedom [22 J. If the time scale of the density (or pressure) 
perturbation is similar to or less than the time scale of these internal relaxation times, 
the temperature response of the system will lag behind that of the perturbation. This 
will prevent the internal degrees freedom from taking up all the heat and will result in a 
decrease in the effective heat capacity]^ This decrease in the effective heat capacity results 
in hysteresis and in dissipation of heat. 

In 1962, Fixman applied the basic ideas of Herzfeld and Rice to describe the viscosity of 
critical mixtures [23 1. He was motivated by the intimate relation between viscosity and at- 
tenuation. Critical mixtures of fluids display a second-order transition which is indicated by 
a critical slow-down of the relaxation rates of the order-parameters. In contrast to Herzfeld 
and Rice, Fixman did not limit his attention to the rates of translational and internal degrees 
of freedom but rather considered a continuum of long-wavelength fluctuations in the order 



"Note that the effective heat capacity will be referred to as the dynamic heat capacity. 



parameter. With this change of perspective he made the connection between the "transfer 
rates" and relaxation rates of order-parameters in viscous systems. The slow-down during 
a transition means large changes in the dynamic heat capacity of the system and thereby in 
the speed of sound. 

Following the argument of Fixman, the slowing-down of the characteristic relaxation rate 
during the lipid melting transition will cause hysteresis and dissipation of heat. Even in 
the absence of critical phenomena, internal friction and heat conduction as introduced by 
Stokes Il24l and Kirchhoff ||25il . respectively, can cause hysteresis and dissipation. However, 
within cooperative transitions these are secondary effects and we will disregard them for 
low frequencies. 

4 System Response To Adiabatic Pressure Perturbations 

Sound is the propagation of a pressure wave that is followed by a temperature wave as a 
consequence of its adiabatic nature. Thermodynamic ally, changes in pressure (dP) and 
temperature (dT) couple to a change in the heat (dQ) of the system: 

where Cp = {dQ/dT)p is the heat capacity at constant pressure and 

Due to a well-known Maxwell relation, {dS/dp)T = —{dV/dT)p, we obtain 

dpJr WJp 



-T 



dS\ fdV 



dTlpKdS,^ 



-cp('f) (15) 



p 



Another Maxwell relation, {dV/dS)p = {dT/dp)s, allows us to write Eq. ( 15 1 as 



dQ\ fdT 



Constant entropy implies that no heat is dissipated into the environment but only moved 
between different degrees of freedom within a closed system. At transitions, the Clausius- 



Clapeyron relatioqjcan be used: 

'dp\ AH 



dT J g TAV 



(17) 



where AH and AV are the enthalpy (or excess heat) and volume changes (excess volume) 
associated with the transition [11|. Note that these are constant system properties for a 
given transition that can be determined experimentally. 

The change in heat (Eq. ( [T3| )) can now be written as: 

dQ = Cp{T,p) (dT - (^) dp] . (18) 



It is clear from Eq. ( 1 8 1 that the heat capacity acts as a transfer function that couples 



adiabatic changes in pressure to changes in heat. 



Eq. ( 1 8 1 governs the equilibrium properties of the thermodynamical system. How- 
ever, here we consider the propagation of sound, which is a non-equilibrium process. 
The theory of sound considers the limit of small changes in pressure and temperature for 
which close-to-equilibrium dynamics can be assumed. This implies linear relations between 
perturbations and system responses. For this reason it is also called 'linear response theory'. 

In any real system transfer rates are finite and changes happen in finite time. Thus, 
the changes in pressure and temperature can be represented as rates, and Eq. ( [T8] ) can be 

rewritten as: 

/TAV\ 



SQ = JdQ = Jcp{t) t-(^^jp 



dt, (19) 



where T = dT /dt, and p = dp/dt are rates. Note that T = TequUibrium^ which holds if 
absolute changes in temperature upon pressure changes are very small. 

If changes in pressure or temperature happen faster than the transfer rate (or relax- 
ation rate), the energy transferred during this change will be only a part of the amount 



otherwise transferred. Considering Eq. ( 18 1, the finite transfer rate will lower the effective 
transfer function, in this case the heat capacity. This means that also the heat capacity 
must contain a relaxation term, (1 — "^cp), with < ^cp < 1- This function describes the 
equilibration of the system. As the system approaches equilibration, (1 — "^cp) approaches 
unity. Below, we will assume that the function ^c^ is an exponentially decaying function 



of time. Eq. ( 19 1 must then be written as a convolution: 

dQ{t) = j [cp(oo) + Acp (1 - m,^{t - t'))] (f{t') - ^p) clt', (20) 



^^The use of the Clausius-Clapeyron relation can be justified by the weak first-order nature of the fipid 



melting transition [261. 



where 6Q{t) is the change in heat, Cp{oo) is the part of the heat capacity that relaxes 
more rapidly than the changes in pressure and temperature considered. In the lipid bilayer 
system Cp(oo) is the heat capacity contribution from lipid chains, which we consider as a 
background contribution. Acp is that part of the heat capacity which relaxes on time scales 
of a similar order or longer than the perturbation time scale. In the lipid membrane system 
this is the excess heat capacity. In Eq. ( [20| ) it has been assumed that the mechanisms of 
relaxation are the same for pressure and temperature. This assumption has been justified 
experimentally and numerically in the literature lfT0l[T4ll27ll28l . 



After partial integration of Eq. pO) ), subsequent Fourier transformation and the use of 
the convolution theorem, Eq. ([20|) can be transformed into (see Appendix): 



/ TAV \ 

6Q = cpico) (t{u) - ^^P{^)] ■ (21) 



T{uj) and p{(jj) can be regarded as periodic variations of temperature and pressure, respec- 
tively. We have also now introduced the frequency dependent heat capacity, 

Cp{u) = cp(oo) - Acp / e-''^*i!c,{t)dt . (22) 

JO 

From Eq. (22i the frequency dependent transfer function (dynamic heat capacity Jj can be 
found, giving a full description of how a lipid bilayer responds to adiabatic pressure per- 
turbations. Both cp (oo) and Acp are experimentally available using differential scanning 
calorimetry (DSC). The only unknown is the relaxation function, ^Cp- 

4.1 Relaxation Function 

The relaxation function of the heat capacity is related to the rate of energy transfer from 
the membrane to the environment. The fluctuation-dissipation theorem ensures that the 
rate of energy transfer is equivalent to the relaxation behavior of energy fluctuations. Since 
the heat capacity is a measure of enthalpy fluctuations, the relaxation function of the heat 
capacity must be the relaxation function of the enthalpy fluctuations ifTTl . 

The relaxation behavior of the fluctuations of enthalpy in pure lipid vesicles has been 
considered theoretically, numerically and experimentally showing that the relaxation of 
enthalpy is well described by a single exponential function ifTOl ITSll : 

[H - {H)){t) = {H- {H)){0) ■ exp (-*-] , (23) 



'^It is important note tiie difference between the dynamic lieat capacity (frequency dependent) and the nor- 
mally known equilibrium heat capacity. The equilibrium heat capacity is a constant system property whereas 
the dynamic heat capacity is an effective heat capacity that can be less than or equal to the equilibrium heat 
capacity as a consequence of the finite transfer rates in real systems. 



where {H — {H)){0) serves only as a proportionality constant and r is the relaxation time. 
For various pure lipid membranes close to melting transitions, it was further found that 
relaxation times are proportional to the excess heat capacity, 



rp2 

T 



T=—Acp, (24) 



where L is a phenomenological coefficient. For LUV of DPPC L = 13.9-10^ J- K/{s-mol) 



4.2 Response Function 

Using the relaxation function of the enthalpy fluctuation as the relaxation function of the 
dynamic heat capacity, 

*cp=expf-^j, (25) 



Eq. ( 22 1 can be solved and the dynamic heat capacity can be determined. 



Cp{uj) = Cp(oo) - Acp / e *'^M — je ^dt 

= c,(oo) + Ac,(^i^) (26) 

Note that the above derivations can be carried out with lateral pressure instead of pressure; 
the choice of using pressure is entirely for notational convenience. 

5 Adiabatic Compressibility 

In estimating the speed of sound in the plane of a lipid membrane during the melting tran- 
sition, the response of the membrane to sound (the dynamic heat capacity) must be related 
to the lateral adiabatic compressibility. The adiabatic lateral compressibility is defined as 

where 11 is the lateral pressure. The adiabatic lateral compressibility can be rewritten in the 
following form ll29l : 

''s=^T-j^[Qf)^, (28) 



where 



''^=-K^).=''^'"'^¥^='' <'" 
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is the isothermal lateral compressibility, k^{cc) is that part of the isothermal lateral com- 
pressibility that relaxes faster than changes in the pressure and temperature considered, and 
^sys em -^ ^^ ^^^^ capacity of the total thermodynamical system, i.e., the lipid membrane 
plus the accessible surrounding aqueous medium that serves as a buffer for heat transfer. In 
the last equality, the empirical proportionality /S.A = ^a^H has been used Il2l|27||, with 
7A = 0.89 m^/J for a lipid monolayer of DPPC. 

In the literature on attenuation and dissipation of sound in critical media a different form of 



Eq. (28 1 is often used to relate the dynamic heat capacity and the adiabatic compressibility, 
using the dynamic heat capacity as the heat capacity of the total system ifTTl |30|| . This 
can be done in a straight forward manner by employing the Pippard-Buckingham-Fairbank 
relations (PBFR) OTllBlJI . The main difference between this approach and the one adopted 
here is that their compressible medium is three-dimensional, and the system heat capacity 
is that of this medium. In contrast, the lipid membrane system is a pseudo two-dimensional 
(the bilayer) embedded in a three-dimensional aqueous medium that serves as a heat reser- 
voir (see Fig. [2]|. Therefore, the aqueous medium contributes significantly to the features 
of the membrane in a frequency dependent manner. 






' \ ' » 



» .^' L 







Figure 2: Visualization of temperature wave in tlie plane of a lipid bilayer. The coloring indi- 
cates heat penetrating into the surrounding water. 

Imagine a standing temperature wave in the bilayer. The transfer of heat from the wave 
to the surrounding water will be time dependent, see Fig. Q for visualization. In the limit 
of cj — )• 0, the amount of water (heat reservoir) participating will effectively go to infinity. 
In the other extreme, (u — )• oo), no heat will be transferred to the surrounding heat reservoir. 
Evidently, the heat capacity of the total system is frequency dependent: 



Cp 



^system ^^^ ^ ^Hpid ^ ^reservoir ^^^ (3^^ 

where Cp"" = Acp + cp (00) is the complete heat capacity (in equilibrium) of the lipid 

11 



membrane and Cp*^''"°*'~(a;) is the heat capacity of the participating heat reservoir. In this 
approach it is the size of the contributing heat reservoir that is frequency dependent. 



Using the proportionality relation A A = ja'^H in Eq. (28 1 and assuming that {dA/dT)ii 
in the chain melting transition region is completely dominated by the transition associated 
change in area, the following approximation can be made [14]: 



A A g*S'''*em 

^ ^(^).'^[Ac,-'pi). (31) 

The parenthesis has the units of a heat capacity and is frequency dependent through the 
frequency dependence of the size of the associated heat reservoir. We pose as an ansatz 
here that this parenthesis is the effective heat capacity of the lipid membrane in a finite 
adiabatically isolated heat reservoir, which is equivalent to the dynamic heat capacity of the 
lipid membrane following the above argument. 

AcpH = Acp-^^g^. (32) 

Cp 

Numerical justification of this ansatz will be published at a later point. 

Using this ansatz, the dynamic heat capacity can be related directly to the adiabatic 



lateral compressibility through Eq. (31 1: 



4 = 4(oo) + ^ • Acp{lo), (33) 

where the Acp{uj) is the dynamic heat capacity without background. In this equation, we 
use the compressibility of the lipid bilayer, which is half of that of the lipid monolayer. 

6 Results - The Speed of Sound 

The goal is to estimate the speed of sound and its frequency dependence in the plane of a 



lipid membrane. From the estimated dynamic heat capacity Eq. (26 1, the adiabatic lateral 



compressibility can found using the proposed relation (Eq. ( 33 1). The lateral speed of sound 
can then be estimated using Eq. (|2]) as: 



c^ 



KJp"^ 
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where k^ is a function of the frequency, w. The effective speed of sound is given by 
(Eq. ^y. 



u 



r,A\2 



Using the previous two equations, one can show that 



n^{u) = ip^) 



A\-l 



Re{KJ) + \ni 



(34) 



Eq. (34 1, the effective speed of sound squared takes the analytic form: 



Inserting the estimated adiabatic lateral compressibility from Eq. ( [33] ) and Eq. (26l) into 

== , (35) 



u'{oj) 



with the notation 



and 



cl + 4 (l + (a.r)2) + 



1_ , 1 

"1 



■2 + q (l + (cJr)2) 



+ 



J^ bJT 

4 (l + ('^r)2) 



cj = {p^4{oo)) 



cl{uj) 



/^Ac,(^) 



(36) 



(37) 



Here, ci is the lateral speed of sound of the membrane outside of the transition, and C2 is 
the component of the lateral speed of sound related to the Upid melting transition. 



All variables in Eq. (35 1 can be found from the excess heat capacity of the lipid melt- 
ing transition and the fluid fractionjj which can be obtained using differential scanning 
calorimetry (DSC). The area, the lateral density and the background isothermal compress- 
ibility are all directly related to the fluid fraction |28|. The relaxation time can be estimated 



from its phenomenological proportionality relation to the excess heat capacity, Eq. (24 1. 



The proposed analytic expression for the effective speed of sound (Eq. ( 35 1) is shown in 
Fig. ([3]), where the excess heat capacity and the fluid fraction are taken from Monte Carlo 
simulations of the lipid melting transition in LUV of DPPC. The simulation have been 
carried out in a manner similar to that described in 1331 . 

The frequency dependence of the speed of sound is described by the function, /{ior) 
with < f{ujT) < 1, defined by 



U^{uJt) =ul + (U^ - ul) f{ujT) 



(38) 



where mq = u{ojt — > 0) and Uoo = u{ojt 



ul 



cxd). From Eq. (35 I we see that 



'1 



The fluid fraction is the fraction of a considered lipid system that is in the fluid phase. 
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Figure 3: Left: The effective lateral speed of sound squared as a function of density at different 
angular frequencies along with the limiting case: u; — >^ 0. Right: The generic function, /(o^r), 
that takes the effective lateral speed of sound squared, at a given lateral density, from the low 
frequency Umit (/(wT — >^ 0) = 0) to the high frequency limit (/(wx —^ (x>) — 1). 



Ur, 



cl 



(39) 



See Fig. ^ (right). The generic function / was chosen to be a function of the dimensional 
less quantity ujt rather than u in order to render it independent of the lateral density. 

6.1 Dispersion Relation 

In the Soliton model described by Eq. Q, dispersion was assumed to be small and in- 
dependent of the lateral density due to the lack of detailed information of the frequency 
dependence of the speed of sound as a function of density. Using the considerations of the 
previous paragraphs, we can now estimate the dispersion in lipid membranes. In the Soliton 
model the extent of dispersion is described by the parameter, h. Assuming that dispersion 
is small, h can be related to the lateral speed of sound as 



Cn + 



huj' 



(40) 
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where Cq ^ u^{oj — )■ 0) = (1/c^ + l/c^) ^. Eq. (WOl corresponds to a Taylor expansion of 
the lateral speed of sound squared to second orderjj Expanding Eq. (35 1 to second order, 



Cq + Cq 



4 3cf + 4ci 



44(4 + 4 



2 2 

-W T , 



(41) 



we see that the dispersion parameter has the form: 






3cf + 4c| 



44(4 + 4) 



(42) 



Using the excess heat capacity and the fluid fraction for large unilamellar vesicles of DPPC 
as used in Fig. Q, we can estimate the density dependence of the dispersion parameter 
h(pA) as shown in Fig. (Q. 
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Figure 4: The dispersion parameter, h, as a function of lateral density for LUV of DPPC, based 
on the proposed expression for the lateral speed of sound. 

Here, the density of the fluid phase is approximately 4 • 10~^ g/m^, the maximum of the 
dispersion parameter corresponds to the chain melting transition maximum, and the density 
of the gel phase is 5 • 10^'^ g/m^. It is clear that the dispersion parameter is strongly depen- 
dent on the lateral density of the membrane. 

The density dependent dispersion parameter, h(p^), will finally enter the differential equa- 
tion (Eq. Q) for the propagating nerve pulse. In the original treatment, h was considered 
an adjustable constant that determined the time scale of a solitary pulse in nerve axons. In 
the present extension, h(p^) is fully determined by the cooperative nature of the membrane 
system and does not contain adjustable parameters. Preliminary calculations indicate that 
this dispersion parameter will yield a natural time scale for the propagating soliton in nerve 
axons. 



*The first order term is zero since the speed of sound squared is symmetric around a; = 0. 
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7 Discussion 

The response of lipid membranes to adiabatic periodic pressure perturbations (sound) is 
closely related to the relaxation behavior of the system ll22l l2?il . Using thermodynam- 
ics and linear response theory, we have described the response of the lipid membrane to 
a perturbation with the assumption that the relaxation function has a simple exponential 
dependence on time. We obtain a form for the dynamic heat capacity which can be under- 
stood as the effective heat capacity when the lipid membrane is subject to periodic adiabatic 
pressure perturbations. The dynamic heat capacity was then related to the adiabatic lateral 
compressibility using the idea that the size of the associated water reservoir is frequency 
dependent |[T4l . The adiabatic lateral compressibility was then used to obtain an expression 
for the effective speed of sound as a function of frequency. 

The major assumption in our approach concerns the nature of the relaxation function. 
We have previously studied the relaxation behavior of the lipid membrane in the vicinity 
of the melting transition at low frequenciesjJThis means that the lipid melting transition is 
assume to be non-critical. The single exponential relaxation behavior should, however, only 
be considered as a low frequency approximation. In a number of ultrasonic experiments 
it has been shown that a single exponential is insufficient to describe the dynamics of the 
cooperative processes involved in lipid melting in the ultrasonic regime |[9l lT4] - [T6l . In these 
ultrasonic experiments, some phase transition phenomena are even apparent in the MHz 
regime. Single exponential relaxation behavior, and thereby the validity of the estimated 
speed of sound, is thus limited to frequencies comparable to the relaxation rate or lower. 

Van Osdol et al. ||T9l have made adiabatic pressure perturbation experiments on unil- 
aminar and multilaminar vesicles of DPPC. They studied relaxation behavior of the lipid 
membrane by measuring the frequency dependence of the effective heat capacity and the 
compressibility as a function of frequency. Although the data available for unilamellar 
vesicles is very limited and has large errors, it can still serve to illustrate qualitative ten- 
dencies of the effective heat capacity, see Fig. Q, that are similar to the theoretical results 
reported here. The effective frequency dependence of the speed of sound shown in Fig. Q 
is dominated by the cooperative properties of the lipid melting transition of DPPC. In this 
model system the relaxation time during the transition is as slow as seconds. In biological 
membranes such as membranes of nerves, realistic characteristic relaxation times can be 
assumed to be of the order of 1 — 100 ms. This change in relaxation times between the 
model system and biological membrane expands the upper limit of the frequency range for 
which our approach is likely to be valid from several Hz to several kHz regime, assuming 
that the general behavior of pure lipid and biological membranes is otherwise similar. Since 
the duration of a nerve pulse is roughly 1 ms, the relevant frequency components contained 
in a nerve pulse can be estimated to be 1 kHz or less. The relevant frequency range for 



'The time resolution of experiments from our group is 0.3 s ~ 3.3 Hz. Relaxation profiles on longer time 
scales are well approximated by a single exponential decav l 1011 181 . 
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Figure 5: Left: The calculated dynamic heat capacity for LUV of DPPC at different frequen- 
cies. Right: The effective heat capacity profiles for LUV of DPPC at different frequencies, 
measured by Van Osdol et al. lHU. The measured effective heat capacities have not been cor- 
rected for contributions from the experimental setup, and a direct comparison is therefore not 
possible. The theoretical dynamic heat capacity shows the same quaUtative features as the 
measurements — a dramatic decrease in the height of the excess heat capacity with increas- 
ing frequency and a relatively constant width. The difference in frequency scales seen in the 
two panels is due to an estimated difference of more than a factor of 10 in the characteristic 
relaxation time. Frequencies are given in units of Hz = {1/2tt) rad/s . 



nerve pulses is thus covered by our proposed expression for the effective speed of sound. 
The present results may thus provide useful insights regarding sound propagation in an 
otherwise inaccessible regime and can extend our understanding of the nature of nerve 
signals. 

In future studies, the hnear response theory described in this chapter will help to define 
an intrinsic length scale of the electromechanical soliton proposed by us as an alternative 
description for the nervous impulse. 

A Derivation of the dynamic heat capacity using the convolution 
theorem 



The purpose of this appendix is to provide additional details in the derivation of the fre- 



quency dependent heat capacity given in Eq. (26 1 starting from Eq. (20 1. The change in 
heat is a convolution of the applied perturbation and the relaxation of the transfer function 
— the effective heat capacity. The perturbation is well defined at all times and can safely 
be assumed to be zero for t — )• — oo. The relaxation function is only defined from [0, c«], 
where t = is the time at which the system starts to equilibrate. The relaxation function, 
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\I', is chosen such that ^(i — ;• 0) = 1 and ^{t — )• oo) = 0. To accommodate the chosen 
form of the relaxation function the convolution can be written as following: 



6Q{t) 



— oo 

t 



{cp{oo) + Acj,{l--^{t-t')))(f{t') 
g{t-t')f{t')dt\ 



TqAF 



p{t') dt' (43) 



(44) 



where g{t — t') is the transfer function and f{t') is the perturbation. Note that 
f{t) = df{t)/dt, Cp{oo) is the component of the heat capacity not associated with 
the melting transition, and Tq is the equilibrium temperature. 



Integration by parts allows us rewrite Eq. (44 1 to the following form: 



5Q{t) 



it 



git') / f{t')dt' 



-oo ^— oo 



f{t")dt"]g{t-t')dt' 



(45) 



The first term in Eq. ( 45 1 takes the form: 



it') / fit')dt' 



/M< 



[9{t')f{t')] 
g{t)f{t)-g{-oc)f{-oc) 



(46) 



where 



fit') = {T{t') - To) 



ToAV 
-^j^{p{t)-Po) 



Assuming that the system is in equilibrium as t — ;■ — oo and f{t' — )■ — oo) = 0, simplifies 

(47) 



Eq. (46l. 



g{t)f{t) - g{-oc)f{-oo) = Cp(oo)/(t). 
The second term in Eq. ( |45| ) can be rewritten by changing the variable to t" = t — t'. 

I' (I f{t')dt'^ git - t')dt' = - j^ fit - t")g{t")dt" , (48) 

where the integration limits have been changed accordingly. 

Since we are interested in sinusoidal perturbations, we consider the Fourier transform 



of Eq. (|43]) and find: 

6Qiuj) -- 



6Q{t)e-''^^dt 

— oo 

oo / POO 

Cp{oo)f{t) + fit- t")g{t")dt" 1 e-"^'dt. 



'iujt , 



(49) 
(50) 



18 



The Fourier transform of the first term in Eq. (50 1 can be carried out without complications: 

POO 

cp(oo) / f{t)e-'^'dt = Cpioo)fiu;). (51) 



The second term of Eq. ( [50] ) can be rewritten as follows: 

9{t") I f{t - t") 



/ f{t-t")g{t")e-'-'dt"dt= / 

-OO Jo Jo 



"\-'^Utdt" 



(52) 



Changing variables again, t' = t — t", the Fourier transform of the second term in Eq. ( 50 1 
can be split into two terms: 



9{t") / fit-t")e-'^'dt 



/•OO /"OO 

t"dt = / git") / f{t')e-'^^''+'"^dt'dt" 

Jo J-oo 







5(t")e"*"* dt" / /(f)e^*"'* dt' 



'iujt • 



g(t)dt. 



(53) 



This is known as the Convolution theorem. From Eq. (53 1 and Eq. ( [5T] ), Eq. (49 1 can be 
written as: 



5Q{uj) = I Cp(oo) + 



'iuit ■ 



where 



g{t)dt) f{oj), 



fiu;) = f{uj) - ^^^P{^) and g{t) = -Acp^(t) . 



(54) 



The Fourier transform of Eq. (|44]) takes the final form: 



6Q{u 



Cp (oo) - Acp 



Cp{u) f{u) 



'iiut 





TqAV 
AH 



4/{t)dt) f{uj 



TqAV 
AH 



p{u}) 



p{uj] 



(55) 
(56) 



Using ^(t) = exp(— t/r), the dynamic heat capacity, Cp{uj), is found to be 

Cpiio) = Cp(oo) + =^ / e-"^'e-'/^dt 



T Jo 



Cp(oo) + Acp 



which has the form of a Debye relaxation term. 



1 — iu)T 

1 + (6^t)V ' 



(57) 
(58) 
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